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Abstract. We show that shearfree perfect fluids obeying an equation of state 
O ! p = (7— are non- rotating or non-expanding under the assumption that the spatial 

divergence of the magnetic part of the Weyl tensor is zero. 
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1. Introduction 



The shear-free fluid conjecture claims that for any general relativistic perfect fluid, in 
which the energy density and the pressure p satisfy a barotropic equation of state 
V = P{^) with p + /i 7^ Cup, necessarily the expansion 9 or the vorticity u vanishes. The 
conjecture has been demonstrated in a number of particular cases: constant p ('dust' 
with a cosmological constant) pill El EE], spatial homogeneity [IT, fTT]. dp/dfi = 1/3 [20] . 
dp/dfi = —1/3 or 1/9 [11 [21], aligned vorticity and acceleration [22], vanishing magnetic 
part of the Weyl tensor H = [7], 9 = #0) [T2J, 9 = 9{u) [H], Petrov types N [3] 
and III [1], models in which there is a conformal Killing vector parallel to the fluid 
flow [6J, purely magnetic perfect fluids [8]. It should be noted that this conjecture 
is not true in Newtonian theory and hence if true in general relativity, then it would 
highlight essential differences, like that of a well-defined universal time, between the 
two theories [17]. Although in most of the above mentioned results a tetrad formalism 
was used, recently it has been shown that a covariant approach can also be successfully 
employed with considerable potential for future analysis [17] [T8l 13]. 

Motivated by the continued absence of a general proof of the conjecture (or lack 
of counter-example), we have decided, as an important step forward, to investigate a 
direct generalisation of Collins' 1984 result that uo9 = holds when the magnetic part H 
of the Weyl curvature (with respect to the fluid congruence) vanishes. In our analysis, 

I otherwise the fluid velocity is not uniquely determined by the geometry and shearfree, rotating and 
expanding examples are known 15j 
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we shall make the weaker assumption that H is solenoidal, that is, that its spatial 
divergence (see below) vanishes. Interestingly, it has been recently shown [T9J [2] that 
the assumption of third order restrictions, such as divH = and/or divE = 0, leads to 
physically interesting families of perfect fluid solutions. Also, in a classification attempt 
of these fluids the shearfree sub-family would appear to be a natural first candidate for 
further investigation. Although a large part of our analysis holds for a general barotropic 
equation of state p = p{fj) (such as the positive definiteness of the matrix and the ensuing 
existence of a Killing vector parallel to the vorticity), at a certain point, we shall need to 
simplify matters by assuming a 7-law equation of state (with a possible non-vanishing 
cosmological term "effectively" present by allowing p = (7 — + constant). In the 
general barotropic case with D b H ab = 0, technical difficulties of a very different nature 
are present and a number of extra subcases remain to be investigated (work in progress). 



2. Tetrad choice and relevant equations 

We consider shear-free perfect fluid solutions of the Einstein field equations 

Rab - ^Rgab = fiU a U b + ptl ab , (1) 

where u is the future-pointing (timelike) unit tangent vector to the flow, /i and p are 
the energy density and pressure of the fluid and h ab = g a b + u a u b is the projection tensor 
into the rest space of the observers with 4-velocity u. The vanishing of the shear can 
be expressed by 

U a -b = \oKb + u ab - ii a u b (2) 

where 9 is the (rate of) volume expansion, u the acceleration and uj is the vorticity. 

Throughout, we will assume familiarity with the notations and conventions of 
the orthonormal tetrad formalism of [14]. Crucial in the successful investigation of 
any problem with tetrad formalisms, is the choice of the tetrad alignment as it can 
dramatically alter the appearance and complexity of the resulting complete set of 
equations. We have chosen our tetrad as follows. We begin by aligning e and e 3 
with u and uj, respectively, such that u> = ue^ ^ 0. The relevant variables then become 
/i, p, u, 6, u a and the quantities Q a , (which determine the rotation of the triad e a with 
respect to a set of Fermi propagated axes), together with the quantities n a p and a a . 
Latin indices will be spacetime indices while Greek indices will take the values 1, 2 and 

3. Further, it is advantageous to replace n a p and a a with the new variables q a ,r a and 
n a , defined by 

n a+ i a -i = (r a + q a )/2, a a = (r a - q a )/2, n aa = n a+1 + n a -i, (3) 

expressions which have to be read modulo 3 (for example a = 3 gives n 12 = (<Z3 + r 3 )/2.) 
The relation between these quantities and the Ricci rotation coefficients can be deduced 
from the commutators fl65l) . We will also use, as extra (extension) variables, the 
components of the spatial gradient of the expansion 

z a = dj (4) 
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and the (3+1) covariant divergence of the acceleration, 

j = u a a = d a ii a + u a ii a - ii a (r a - q a ). (5) 

As neither the Jacobi identities nor the field equations contain expressions for the 
evolution of the Q a , it is good practice to choose the triad e a such that fl + uj = 0. The 
fact that this is always possible follows from the results of the action of the commutators 
[83, di] and [83 , 82} on p. One finds that Qi — Q2 — after which a rotation in the 
(12)-plane can be chosen such that fl 3 + to = 0. Herewith, the tetrad is determined up 
to rotations in the (12)-plane over an angle a satisfying <%« = 0. Noticing that the 
evolution equations for the quantities n u — n 2 2 and n 12 are identical (see [H]), allows 
one to further fix the tetrad by making either ran — 1122 = or = 0. Henceforth our 
choice will be 

nu = n 22 = n. (6) 

We express the vanishing of the spatial divergence of the magnetic part of the Weyl 
tensor, H ab = C ac bdU c u d , by means of the Bianchi identity [T3j . 

(divH) a = D b H ab = 3cu b E ab + (jm + p)u a . (7) 

Here the spatial derivative operator is defined by 

D a S c '" d e ...f = h a b h c p ■ ■ ■ h d q h e r ■ ■ ■ hf S VbS p "' 9 r ...s- 

This shows that when D b H a i > = 0, the vorticity is an eigenvector of the electric part 
of the Weyl tensor, E ab = C ac bdU c u d , with eigenvalue — (/i + p)/3. With our choice of 
tetrad this implies 

E 13 = E 23 = 0, (8) 
£33 = -£±* (9) 

In order to relate the components of E ab with the spatial gradient of the acceleration, 
we will make use of the Ricci identity [13] 

E a b = D {a il b) - LO {a UJ b) + U( a U b ), (10) 

where S( a b) stands for the spatially projected and trace-free part of S ab . 

The basic equations of the formalism are now the Einstein field equations and the 
Jacobi equations, which we present, using the simplifications above, in the appendix. 
Henceforth we will assume T^lpllCEHllS], 7 ^ f £0], 7 ^ § [2], 7 ^ f gl] and, of 
course 7 7^ (cf. introduction). 

We shall argue by contradiction in order to establish that a shearfree fluid under 
the given conditions satisfies uj6 = 0. The important step is contained in section 3, 
where we prove first that, if oj6 7^ 0, a Killing vector exists parallel to the vorticity. This 
results in a great simplification of the governing equations and, although this subcase 
still contains all the intricacies of the 'general' problem, it allows us to complete the 
proof in section 4. 
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3. Proof of the existence of a Killing vector parallel to lj for ujO ^ 



First notice that the equations (173 II 76 j) immediately lead to evolution equations for the 
variables r a and q a , 

3<9 r a = —9u a — z a — 9r a (11) 
3d q a = z a + 9(ii a + q a ) (12) 

while ( 1661) and the (Oa) field equations ( J7H1I801) give us the spatial derivatives of uj 
2 

diu = -z 2 - + 2ui) (13) 
2 

<9 2 cj = —-2! + uo(r 2 - 2ii 2 ) (14) 

d 3 u = u(ii 3 + r 3 - g 3 ) (15) 

together with the algebraic restriction 

2 

™33 = ^-*3- (16) 
3cj 

The evolution equation for n follows from (l74j) : 

Q 

d n = --n. (17) 

Acting with the commutators [do, d a ] and [d±, d 2 ] on the pressure and using ( |72l) together 
with the conservation laws 

doV = -9(Li+p) (18) 

<9 Q p = ~(n + p)u a (19) 

leads to a first set of evolution equations for the acceleration and vorticity: 

4 

d u a = (7 - l)z Q + (7 - v)9u a , 

5 3 (20) 

8 uj = (7 - -)w0. 

The spatial derivatives of the acceleration can be obtained from ffTU]) . using (jSJ): 



= 


1 2 1 2 n 

--^ + -J - «2<?2 + «3?"3 - «i + -Ell 


(21) 


d 2 u 2 = 


1 2 1 2 n 

-gW + -J — %g 3 + u x n - m 2 + £22 


(22) 


d 3 u 3 = 


1 2 1 2 n 

-gW + -j — wigi + w 2 r 2 - % + £ 33 


(23) 


diii 2 = 


u;0(l - 7) + q 2 ii\ + ^n 33 w 3 - Wi« 2 + #12 


(24) 


d 2 iii = 


-io9{\ - 7) - riw 2 - 2 n 33«3 - ^1«2 + #12 


(25) 


dtu 3 = 


-7*3Wi - 7^2^33 - UiU 3 + Ei 3 


(26) 


d 2 u 3 = 


7^1^33 + <?3« 2 - M 2 «3 + #23 


(27) 
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d 3 ux = -^u 2 n 33 + nu 2 + q\U 3 - U1U3 + E 13 (28) 

1 . 
d 3 u 2 = 2^1^33 - nui - r 2 u 3 - u 2 u 3 + E 23 (29) 

Next we act with the [do, d a ] commutators on u and 9 and use the propagation of (fTE!) 
along u in order to obtain expressions for the evolution of z a and the spatial gradient 
of j: 

9t — 10 

d oZl = 0( 7 - 2)2! - u (z 2 + 9u 2 ) 

9 ^ 10 ll ,(,,+9uA (30) 



d z 2 = 0(7 - 2)^2 H u (zi + 9u\ 

doz 3 = 0( 7 - 2)z 3 



and 



<?lj = 0*1 (7 - 1) - Z 2 U - UxJ - - Mi// — U]J) 

+ -m 1 (18cj 2 + # 2 ) -u 2 — 0u; + 4u; 2 gi (31) 

«, . fl , -14 + 277 ■ ■ -l + 2 7 . 7-1. 
o 2 j = ^2(7 - 1) + z x u - U 2 J - - — U 2 fM — u 2 p 

I Qry _ 10 

+ -ii 2 (18u 2 + 9 2 ) +«!— 9uo-Auo 2 r 2 (32) 

9 3 j = ^3(7 - 1) - u 3 j - 2^-1) ^ ~ ^~2~^' dP + ^3(-18u; 2 + 9 2 ) 

- 4(r 3 - q 3 )oo 2 . (33) 

Now we may evaluate ^2 a [do, d a ]u a by using (jSJ) and ( 12111231) and this leads to the 
evolution of j as 

d J = (7 - l)^ a + 0(7 - ~)j + (27 - !j)^ a + (4 7 - y 

+ ( 7 -l)(g a -r Q )^ (34) 

Next a long calculation, involving the propagation of the field equations ( 18TH861) along 
u and the use of the [d , d a ] commutators on rp,qp,itp together with the [d a , dp] 
commutators on 9 and lu, allows one to obtain algebraic expressions for the directional 
derivatives d a zp, as well as evolution equations for j and E a p (the latter can equally well 
be obtained by writing out the 'dot E' Bianchi identities [13]). The resulting expressions 
can be simplified by the introduction of the following non-linear combinations of 
kinematic quantities, each having a clear geometric meaning, 

O O / \ 

U = % + w 2 , V = U1Z1 + u 2 z 2 , W = u\z 2 — u 2 Zi, Z = z 1 + z 2 . (35) 

We obtain: 

a) the evolution of j: 

d J = + 0U + z 3 u 3 + 9u 2 3 ) + (7 - 1)(9 7 - 1O)^ 2 + ^-^'0; (36) 
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b) the evolution of E a/3 : 
doEu = — -9E n 



<9o-&22 — — ^9E 22 



^ (ui*i - 2m 2 z 2 + u 3 z 3 - (cu 2 (9 7 - 8) - til + till - iil)6), (38) 



9(3 7 - 2) 



d E 12 = ~9E 12 + \, (u2Zi + uiz 2 + 2miM 2 ^), (39) 

2 

d £i 3 = ~-0E 13 + 7 !,. ^ (^i + «i* 3 + 2miM 3 ^), (40) 



3(3 7 


-2) 


3 7 - 


-4 


3(3 7 


-2) 


3 7 - 


-4 


3(3 7 


-2) 



2 . . . 

d E 23 = --9E 23 + 777 ^- — ^(u 3 Z2 + u 2 z 3 + 2u 2 ii 3 9); (41) 

c) the spatial derivatives of z a : 

2(15 7 ~ 14) . ^ 2 (6 7 -8)+g 2 (6-97) 
= - 3(3 7 - 2) UlZl + 3(^2) Z2 + Z3 

+ 3(3 7 ^-2) (Lj2(2 " 337 + 2?72) " (67 " 8)(2 ^ " "2 " "a) + (6 " 97)^n), (42) 

2 1 07 - 14) . , M 

<7 2 ,2 2 = ^^2-22 H ; s z i H ; ; ^3 

3(37 - 2) 3(37 -2) 

+ 3(3 7 -2) (a;2(2 " 337 + 2?72) + (67 " 8)( " 2 " ^ + ^) + ( 6 " 97 )^ 22 )' ^ 43 ) 

2(15 7 - 14) . , , 

dszs = - 3(37 _ 2) u 3 z 3 + ^— ^ Zl + z 2 



^3(67 


-8)+r 3 (9 7 -6) 




3(3 7 - 2) 


!) + (6 


-9 7 )£ n ), 


^3(67 


-8)-g 3 (9 7 -6) 




3(3 7 - 2) 


:) + (6 


-9 7 )E 22 ), 


?i 2 (67 


-8) -r 2 (6-97) 


3(3 7 - 2) 
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s . rK(56 - 78 7 + 27 7 2 ) + (67 - 8)U + (16 - 12 7 )7j 2 + (6 - 9 7 )£ 33 , (44) 
3(37 - 2) 

d x z 2 = q 2 zi - ^p- — ^-{u 2 zi + u x z 2 ) + \n 33 z 3 - ^2 — ^-u x u 2 9 - 9E 12 
07 — z z 07 — z 

+c(^- J + V-2^ 2 ), (45) 

<9 2 ^i = -riZ2 - ^ — ^(^2^1 + u x z 2 ) - \n 33 z 3 - 2 ^ 3 - 7 - — ^-U\U 2 9 - 9E 12 
07 — z z 07 — z 

--(^-J + V-2- 2 ), (46) 

cfczi = q x z 3 - i(n 33 - 2n)z 2 - ^ 7 _ ^ (u 3 z 1 + uiz 3 ) - 2 ^ 37 _ ^ U\u 3 9 - 9E 13 , (47) 

1 6(7-1) 2(37-4) 

<9 3 ^ 2 = ~r 2 z 3 + o( n 33 - 2n)zi - — ^-(^3^2 + U2Z3) — u 2 u 3 9 - 9E 23 , (48) 

z 07 — z 07 — z 

1 6(7-1) 2(37-4) 

diz 3 = -r 3 zi - -n 33 z 2 - — - - (tt 3 zi + uiz 3 ) — - - ttitt 3 g - 9E 13 , (49) 
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1 6(7-1) 2(37-4) 
d 2 z 3 = q 3 z 2 + -n 33 zi - — — (u 3 z 2 + u 2 z 3 ) 5 ^-u 2 it 3 - 0E 23 . (50) 

L 07 — I 07 — I 

The evolution of the new variables U, V, W and Z is given by 

9of / = !P2_zi) w + 2(7 - 1) v 

9oV ^J^ev-^ w + b _ l)z <5i) 
m = h^ e .u + ?X^v + 2 < 3 ^ 5 W ' 

° 2 2 3 

d Z = 0lu(9j - 10)W + 2(7 - 2)9Z 

The equations given so far describe the behaviour of a general shearfree perfect fluid 
with a 7- law equation of state (and hence can be useful in further investigations of the 
ujQ = conjecture). It is only at this stage that we introduce the condition that the 
magnetic part of the Weyl tensor is divergence-free. Imposing ([HI M) we deduce from 

(7 - 3) (z 3 + 29u 3 )ux + u 3Zl = (52) 

(7 - (*3 + 29u 3 )ii 2 + u 3 z 2 = 0, (53) 

where the factor 7 — | can be divided out (7 = I corresponds to a pure radiation perfect 
fluid, in which case the conjecture uj6 = is known to hold [20]). 

3.1. z 3 + 26u 3 ^ 

First suppose that z 3 + 26u 3 ^ 0. Now propagating equations (I52II53P along u results in 
a linear system in zi, z 2 , 

az\ + (3z 2 = 0, —/3zi + az 2 = (54) 



with 



a = (7 - l)(z 3 + u 3 9f + ( - j + (7 - V - 2u 2 ) ul 



2 J w 3' 



97 — 10 
P = 7 (zs + 0u 3 )uu 3 . 



(55) 



If /? 7^ the system (1541) has only the 0-solution, and consequently (I52II53P would imply 
iii = ii 2 = 0. In this case, the acceleration would be parallel to the vorticity and the 
conjecture uj9 = follows by [22]. When (3 = we have the following possibilities: 

• 7 = y : i n this case the conjecture u6 = has been demonstrated [21] 

• u 3 = 0: (1521 now implies % — ^2 = which takes us back to the 'dust' cases with 
cosmological constant [HI EDI EE] 

• z 3 + 6u 3 = and hence, by ( 1521) . 1*3(21 + Qui) = ii 3 {z 2 + 6u 2 ) = 0. Assuming u 3 ^ 
this implies that the spatial gradient of log# — J vanishes, and the fluid flow is 
irrotational (unless of course \og9-J = 0, in which case the u9 = conjecture 
would follow from pj]), which contradicts our assumption. 
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3.2. z 3 + 29u 3 = 

If z 3 + 29us = then ( l52|) implies that w 3 ,2i = -u 3 z 2 = 0. If -ii 3 7^ propagation of the 
conditions z\ = z 2 = along it gives, using ( 1301) . 

Ui(97-10)^ = m 2 (97-10)^ = (56) 

which again would imply that the vorticity and acceleration are parallel. We therefore 
conclude that ii 3 = and hence also z 3 = 0. By (IBUl) one also has then n 33 = 0. 
Furthermore, with (I26II27I) and (I49II50P the conditions ii 3 = z 3 = yield 

r 3 z x = q 3 z 2 = 0, r 3 Mi = q 3 ii 2 = 0, (57) 

such that, unless acceleration and vorticity are parallel, at least one of q 3 or r 3 = 0. 
Without loss of generality we can suppose q 3 = 0. If r 3 7^ then z\ = iii — 0, after which 
( 130]) would imply (97 — 10)07(2:2 + 6u 2 ) = 0. The spatial gradients of 6 and p would 
then be parallel (with the coefficient of proportionality depending on 9 only), which 
again would imply that the flow is irrotational. We conclude that both q 3 = r 3 = 0. A 
straightforward calculation shows then that there is a Killing vector along e 3 , although 
in the sequel this property will not be explicitly used. 



4. Proof of ojQ = 

From (I23II44I) and ii 3 = z 3 = one obtains 

j + 2lo 2 - 3u iqi + 3ii 2 r 2 + 3E 33 = (58) 

(3 7 - 2)(-3q lZl + 3r 2 z 2 - 36E 33 ) + (3 7 - 4) (26 U + 2V + (9 1 - U)uo 2 6) = (59) 

Propagating (Q along u gives a third algebraic relation among the same variables, which 
allows us to express qi and r 2 as 

qt = ((7 -l)(p + /j l ) + 3(3 7 - 4V 2 ) Qii 2 W- x + (j-p-/i + 2oo 2 )z 2 (3W)- 1 (60) 

r 2 = ((7 - l)(p + fi) + 3(3 7 - 4V 2 ) einW" 1 + (j-p-fi + 2u 2 )z l (3W)- 1 (61) 

and by which we can rewrite (1591) as 

(37 - 2) 2 (p + fj)6 + (37 - 4) (4(97 - 8)9lu 2 + 2V + 28U) = 0. (62) 

Note that W 7^ since propagating W — along u would imply (97 — 10)Uu(zi + 9ui) = 
and hence z\ + 9iii = z 2 + 8u 2 = 0. Applying the same reasoning as in section 13.21 
above would then lead to an irrotational flow. 

We now focus on equations (I31|32p . Applying the [di, d 2 ] commutator on j and 
using (PI El EU E21 M, 12511421 SSI HSl \M EHl ED) we obtain 

2(4 7 -5)l^+(77-10Vj + (37-2)(ll7-8)(p + /i)u;+(4597 2 -10487 + 580V 3 = 0.(63) 

Propagating this relation twice along u and simplifying the result using ( 1591) and ( 1631 ) 
results in a linear system 

67(217 - 20) (37 - 4) (37 - 2) (477 2 - 767 + 30) 
4(5 - 37) (2I7 - 20) (37 - 4) (37 - 2) (477 2 - 767 + 30) 

which clearly implies u; = 0orp + yu = 0. 



p + p, 



(64) 
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We have demonstrated that a shearfree perfect fluid, obeying an equation of state 
p = (7— and satisfying the 'solenoidal' condition divH = 0, is either expansion free or 
vorticity free. Our proof relies heavily on the subcases discussed earlier in the literature 
for various values of 7. In addition, we have presented an interesting framework which 
shows promise for successfully mounting an assault on the shearfree conjecture with the 
assumption of just a gamma law. Presently we are trying to generalise our result to 
include an arbitrary barotropic equation of state p = p(p) with the solenoidal condition. 
Preliminary work shows that almost all relations in the present work can be generalised 
(in particular the fact that u6 ^ would imply the existence of a Killing vector parallel 
to the vorticity). Some technical difficulties still remain to be dealt with in the final 
stage of the proof. This is primarily due to the resulting more complicated nature of 
equation fl63|) and its subsequent derivatives. A proof of the conjecture would be very 
desirable, as the classification of the shearfree case would form a necessary and natural 
first step in the study of purely 'solenoidal' perfect fluids. 



6. Appendix 

Commutator relations, using a a p = 0: 

[do, d x ] = iiido - 6>i<9i + (uj 3 + tt 3 )d 2 - (u 2 + tt 2 )d 3 

[d , d 2 ] = u 2 d - (0J3 + $h)di - 6 2 d 2 + (wi + Vt{)d 3 

[do, d 3 ) = u 3 d + (u 2 + ton)di - 9 3 d 3 - {u x + Qi)d 2 

[di, d 2 ] = -2uj 3 d + q 2 di + r x d 2 + n 33 d 3 

[d 2 , d 3 ] = -2^0 + q 3 d 2 + r 2 d 3 + n n d x 

[d 3 , dx] = -2uj 2 d + r 3 di + qid 3 + n 22 d 2 



Using in addition the simplifications uoi = uo 2 = Vti = Vt 2 = 0, Vt 3 = —uo 3 , n\\ = n 22 = n 
one obtains the: 

Jacobi equations: 

d 3 u = u(u 3 + r 3 - q 3 ) (66) 

din + d 2 r 3 + d 3 q 2 - 72(77 - Qi) - r 3 r 2 + q 3 q 2 = (67) 

d 2 n + d 3 n + d x q 3 - n(r 2 - nq 2 ) - r 3 ri + q 3 q x = (68) 

2 

d\r 2 + d 2 q x + d 3 n 33 - ~6u - r 2 r x + q 2 q x - n 33 (r 3 - q 3 ) = (69) 

d 2 ii 3 - d 3 ii 2 - nui - q 3 ii 2 - r 2 ii 3 = (70) 

d 3 ii! - dxii 3 - nu 2 - r 3 u x - q x u 3 = (71) 

4 

2d u + d\ii 2 - d 2 ii\ - ii\q 2 - ii 2 rx - ii 3 n 33 + -Qu = (72) 
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d (r a - q a ) - -d a 9 + z a + -(r a - q a + 2u a ) = 

3d n + 3d 3 u + n6 — 3oj(u 3 + r 3 — q 3 ) = 
3d n 33 + 3d 3 u + n 33 9 - 3u(ii 3 + r 3 - q 3 ) = 

<9 (r a + g a ) + -0(r a + q a ) = 



Einstein equations: 

«9o£ = -^ 2 + 2^ 2 -i(/i + 3p)+j 
2/3^1 + d 2 uj — a;(r 2 — 2u 2 ) = 
2/3,22 — d\uo — uo(qi + 2ui) = 
2/3z 3 - ujn 33 = 

-<9ir 2 + 9 2 gi - rir 2 - gig 2 - 2r 2 g x - 2r 3 n + n 33 (r 3 + q 3 ) - 2q 3 n = 

-diii 2 - d 2 ui - 2ii l ii 2 + u x q 2 - ii 2 r 1 
-d 2 r 3 + d 3 q 2 + din - <9in 33 - n{r\ - g x ) - 2gin 33 - r 2 r 3 - q 2 q 3 - 2r 3 q 2 = 

-d 2 u 3 - d 3 u 2 + ?ii(n 33 -n) - 2u 2 u 3 + u 2 q 3 - ii 3 r 2 
diq 3 - d 3 r 1 + d 2 n 33 - d 2 n - r 2 (2n 33 - n) - nq 2 - - q x q 3 - 2r x q 3 = 

-<9i« 3 - 9 3 Mi - Mi(2m 3 + r 3 ) + ii 2 (n - n 33 ) + u 3 q x 

-d 1 r 1 + + <9 2 g 2 - d 3 r 3 + q 2 + r\ + -n 2 33 - nn 33 + r\ + q\ - r 2 q 2 - r 3 q 3 
2 8 1 

9° 2 + 3^ ~ 3 ^ ~ ^ ~ dliil ~~ + ^ 3ra ~ ^ 2<?2 
-d 2 r 2 + d 2 q 2 - d 1 r 1 + <9 3 g 3 + q 2 3 + r\ + ^n 2 3 - nn 33 + r 2 + g 2 - r x q\ - r 3 q 3 
2 8 1 

-6» 2 + -u 2 --(2/j,- j) - d 2 ii 2 -u\ + iixn - ii 3 q 3 

-d 3 r 3 + d 3 q 3 + diqi - d 2 r 2 + q\ + r\ - \r? 33 + r 2 + q\ - r x q x - r 2 q 2 = 
2 1 

-9 2 + - (2uj 2 - 2/i + j) - <9 3 m 3 - m 2 + u 2 r 2 - iixqi 
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